arXiv:hep-ex/9706028v2 1 Jul 1997 


D0 Results on W Boson Properties 


Kathleen Streets* 

(for the D0 Collaboration) 

New York University, NY, NY 10003, USA 

(February 7, 2008) 


Abstract 


The D0 experiment collected Ri 15 pb~^ in run lA (1992-1993) and « 
89 pb“^ in run IB (1994-1995) of the Fermilab Tevatron Collider using pp 
collisions at ^/s = 1.8 TeV. Results from analyses of events with W and 
Z bosons are presented for the run IB data samples. From W ev, pv 
and Z —> ee, pp decays, the W and Z production cross sections and the W 
width are determined. Events with W ^ tv decays are used to determine 
the ratio of the electroweak gauge coupling constants as a measure of lepton 
universality. Using W ^ ev and Z ^ ee decays, the W boson mass is 
measured. 


*Presented at the Hadron Collider Physics XII Conference, 
June 5 - 11, 1997, Stony Brook, New York, USA. 
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I. INTRODUCTION 


The D0 experiment collected ~ 15 pb“^ in run lA (1992-1993) and ~ 89 pb“^ in run IB 
(1994-1995) of the Fermilab Tevatron Collider using pp collisions at y/s = 1.8 TeV. Results 
are presented from data collected by the D0 experiment that test the Standard Model 
(SM) of electroweak interactions |^. Measurements of the W and Z boson production cross 
sections, the W decay width, the ratio of the gauge coupling constants, and the W mass are 
presented. 


II. THE D0 DETECTOR 

The D0 detector was designed to study a variety of high transverse momentum (p^) 
physics topics and has been described in detail elsewhere [@. It does not have a central 
magnetic held, making possible a compact, hermetic detector with almost full solid angle 
coverage. The detector has an inner tracking system which measures charged tracks to a 
pseudo-rapidity rj < 3.2, where rj = — In tan | and 9 is the polar angle. The tracking system 
is surrounded by hnely-segmented uranium liquid-argon calorimeters (one central and two 
end-caps). Surrounding the calorimeter is a muon magnetic spectrometer which consists of 
magnetized iron toroids that are situated between the hrst two of three layers of proportional 
drift tubes. 

Electrons and photons were identihed by the shape of their energy deposition in the 
calorimeter and a matching track (for electrons). The energy (E) was measured by the 
calorimeter with a resolution of ~ 15%/(GeV). Neutrinos were not identihed in the 
detector but their transverse momentum was inferred from the missing transverse energy 
in the event: where the sum i extends over all cells in the calorimeter. 

Muons were identihed by a track in the muon chambers matched with a track in the central 
tracking chambers. 


III. IT AND Z PRODUCTION 

Events in which a IT or Z boson is produced are used to measure the cross section times 
branching fraction, the IT width and the ratio of the gauge couplings. In these analyses, the 
IT and Z gauge bosons are identihed through their leptonic decay modes: IT —>• eu, pu, tv 
and Z —>■ ee, pp. These modes have a cleaner signature and are easier to distinguish from 
the background of QCD multijet production than hadronic decay modes. The events with 
decays into e’s and p’s are selected by requiring a high-pr e or p and large $rp for IT’s and 
two high-pT e’s or p’s for Z’s. The hadronic decay of the r is used to to select the W ^ tv 
events. 


A. Production Cross Sections 

The measurement of the product of the cross section and the branching fraction for IT’s 
and Z’s provides a fundamental test of the Standard Model. These measurements have been 
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published for the run lA data sample and the preliminary results are presented here for 
the run IB data sample. 

For the hnal event selection in this analysis, electrons were restricted to a region \r]\ < 1.1 
and 1.5 < \t]\ < 2.5 and muons to a region \r]\ < 1.0. The IF —> ez/ events were selected 
by requiring the transverse energy of the electron Et > 25 GeV and > 25 GeV and 
the Z —>■ ee events were required to have two e’s with Et > 25 GeV. The IF —>■ /rz/ event 
selection required > 20 GeV, and I^t > 20 GeV and the Z —> /i/r selection required 

Pt > 15, 20 GeV for the two /r’s. The transverse mass for IF events and invariant mass for 
Z events in the hnal data samples are shown in Fig. |^. Table | gives the number of events 
observed, the acceptance, the efficiency, the background and the luminosity for these data 
samples. 
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FIG. 1. Transverse and invariant mass distributions for the IF —> ev,p,v and Z —> ee, pp, run 
IB data samples. 
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Electron 

Muon 

# W candidates 

59579 

4472 

Acceptance (%) 

43.4 ± 1.5 

20.1 ±0.7 

ew{%) 

70.0 ± 1.2 

24.7 ± 1.5 

Background(%) 

8.1 ±0.9 

18.6 ±2.0 

Luminosity {pb~^) 

75.9 ± 6.4 

32.0 ±2.7 

^ Z candidates 

5702 

173 

Az{%) 

34.2 ±0.5 

5.7 ±0.5 

ez{%) 

75.9 ± 1.2 

43.2 ±3.0 

Bkg.(%) 

4.8 ±0.5 

8.0 ±2.1 

Lum. {pb~^) 

89.1 ± 7.5 

32.0 ±2.7 


TABLE I. The quantities used to measure the preliminary cross sections for W ev, and 
Z —> ee, for the run IB data sample. 


The preliminary measurements of the cross section times branching fraction [a ■ B) are 
given in table |I| and are shown in Fig. ^ along with the results from CDF [Q. The r results 
shown will be discussed in section [111 Q . Also shown in Fig. |] are comparisons of a ■ B with 
SM predictions |^. The predictions use the CTEQ2M parton distribution functions (pdf) 

i- 



e ^ e n I e 


FIG. 2. Tevatron measurements for the cross sections times branching ratios for W —> eu, /m, rn 
and Z —> ee, pp compared to SM predictions. 


B. W width 

The ratio of the W and Z production cross sections can be used to measure the leptonic 
branching ratio B{W —>■ lu) and extract the W width (Fw)- From the measured width, 
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aw ■ B{W 

—> ev) 

= 

2.38 

± 0.01 

± 0.09 

± 0.20 nb 

aw ■ BiyV 

pj^) 

= 

2.32 

± 0.04 

± 0.16 

± 0.19 nb 

az ■ B{Z - 

-> ee) 

= 

0.235 

± 0.003 

± 0.005 

± 0.020 nb 

az ■ B{Z - 

IJ-IJ) 

— 

0.202 

± 0.016 
(±stat) 

± 0.020 
(±syst) 

± 0.017 nb 
(±lum) 


TABLE II. 

The preliminary cross 

sections for W - 

-> ev, gv and Z — 

> ee, gg. 


a limit may be placed on unexpected decay modes of the W. Many common systematic 
errors, including the luminosity error, cancel in the leptonic branching ratio: 

B{W Iv) _ aw T{W Iv) Fz 
az-B{Z^ll) ~ az T{Z^ll) Tw' 

Using the results above for a ■ B and combining the electron and muon measurements, we 
obtain a preliminary run IB result of i? = 10.32 ± 0.43. The leptonic branching fraction of 
the W may then be calculated, B{W —> lu) = B{Z ll) ■ {az/aw) ■ R = (10.43 ± 0.44)% 
using the measured value of i?, the value oi B{Z ll) from LEP measurements and 
= 3.33 ± 0.03 from the SM prediction [§|. The total width of the W is then obtained 
from this measurement of B{W —> lu) and the value of T{W lu) from SM predictions 0. 
The preliminary run IB measurement is 

Tw = 2.159 ± 0.092 GeV. 

Comparison of the published world average Tw = 2.062±0.059 GeV [Q (does not include 
the run IB measurement) with the SM prediction Tw = 2.077 ± 0.014 GeV p, gives a 95% 
conhdence level upper limit of AT^y < 109 MeV on unexpected (non-SM) decays of the W. 


C. Measurement of the Ratio of the Couplings 

The decay hU —rz/ is studied as a test of lepton universality by measuring the ratio 
of the electroweak coupling constants g^/gY- The r events are obtained from a sample 
in which inelastic collisions were selected by requiring a single interaction signature from 
the Level 0 trigger. The integrated luminosity for the r trigger used in this analysis is 
16.8 ±0.9 pb-h 

To select the W —>■ ru events from the W data sample, the hadronic decay of the r 
is used. These events are identihed by the presence of an isolated, narrow jet. Jets were 
reconstructed using a cone algorithm with radius 0.7 in 77 — 0 space and the width of the 
jet was required to be rmsjet < 0.25. The requirements that E 7 ’(jet)> 25 GeV, {\ri\ < 0.9), 
$rp > 25 GeV and that there be no opposite jet were placed on the data sample. In order 
to separate the events with a jet from a r decay from the large background of QCD jets, 
the prohle distribution of the jets is used. The prohle is dehned as the sum of the highest 
two tower divided by the cluster Et- The prohle distributions from the r sample and 
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the QCD background sample (selected from events with low Ipj, ) are shown in Fig. A 
requirement that the prohle variable be > 0.55 is made to select the hnal r event sample. 
The shaded low-prohle region in Fig. ^ is used to estimate the remaining QCD background. 




FIG. 3. The distribution of the profile quantity for the r candidate sample and for the QCD 
background sample. 


The number of signal events contained in the hnal data sample is listed in table pi 
along with the estimated background contributions. The preliminary value of the cross 
section times branching ratio is a ■ B{W —>■ tu) = 2.38 ± 0.09(sfat) ± 0.10(s|/st) nb where 
the error due to the luminosity has not been included. Comparing this value with the 
measurement of a ■ B{W —> ev) from run lA |^, the ratio of the couplings is determined 
g'^I qY = 1.004 ± 0.019(stat) ± 0.026(s|/st). This measurement shows good agreement with 
e — r universality at high energy. 


IV. W MASS 

The electroweak Standard Model can be specihed by three parameters. These may be 
taken to hea, and Mz, all measured to < 0.01%. At lowest order, the W mass is precisely 
dehned as M\y = A/ sin 9]y with sin^ 9^ = 1 — where 9w is the weak mixing angle 

and A = {j^aj■ The current data are sufficiently precise to require comparison 
to theoretical predictions which include higher order corrections. These corrections have 
contributions due to the running of a —>■ a{M‘^) and to loop diagrams which introduce a 
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Number of Events 

Final Data Sample 

1202 

QCD Background 

106 ± 7 ± 5 

Noise Events 

81 ± 14 

Z ^ TT 

32 ±5 

W 

3± 1 


TABLE III. The quantities used to measure the run IB preliminary cross section for W tv. 


dependence on the sqnare of the top qnark mass, rritop, and the log of the Higgs mass, Mh- A 
precision measnrement of the W mass therefore dehnes the size of the radiative corrections 
in the SM and along with rritop it can constrain Mh- A direct measnrement also serves to 
test the consistency of the SM. 

Previous results from the run lA data sample have been published |]^ and yielded a 
value of Miy = 80.350 ± 0.270 GeV/c^. In the analysis presented here, the preliminary 
measurement of M\y from the run IB data sample is presented, using a calorimeter-based 
measurement. The calorimeter is not calibrated independently to the precision needed and 
therefore the ratio of the W to Z masses was measured and then scaled to the precisely 
known (< 0.01%) Z mass Many systematic errors cancel in this ratio. 

Experimentally, the remnants of the interaction pp (—^ ev) -|- X are detected. Here 

X is due to the recoil (rec) to the W plus the underlying event. The energy of the electron 
and the were measured. The $rp = —prp[rec) —Pt{^) = Pt{W) —prie) and is identihed 
with the neutrino transverse momentum Pt{i^) but differs fromp 7 ’(z/) because of the presence 
of the underlying event. 

Because the longitudinal momentum of the v is not measured, the W invariant 
mass cannot be constructed. Instead the distribution in transverse mass M^iW) = 


^j2pT{e)pT{y) — ‘^Pt{g) ■ Pt{^) is used to obtain the W mass. For Z decays, the energies of 
both electrons are measured and the invariant mass is reconstructed. 

The W ^ ev events were selected by requiring an isolated electron with Et > 25 GeV, 
PriW) < 15 GeV/c and $rp > 25 GeV. The Z —>• ee events were selected by requiring 
two isolated electrons each with Et > 25 GeV, and 70 < Mz <110 GeV/c^. Electrons 
were required to be in the region \p\ < 1.0. There were 28323 W events and 2179 Z 
events in this sample. The electron polar angle was determined from the shower centroid of 
the energy cluster in the electromagnetic (EM) calorimeter and the center-of-gravity of the 
corresponding track. The uncertainty in determining this angle results in an uncertainty of 


±28 MeV/c^ in M^y. 

The mass of the W is determined by a maximum likelihood ht of the measured Mt{W) 
distribution to Monte Garlo (MG) distributions which were generated for 21 different values 
of Myr in 100 MeV steps. This fast MG simulation uses a theoretical calculation for the W 
production and decay and a parameterized model for the detector response. Kinematic cuts 
are placed on the MG quantities as done in the data. All the parameters in the MG are set by 
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FIG. 4. The invariant mass distributions are shown for (a) vr® —> 77 —> e^e“e^e“ events 
(points), (b) J/V’ —> ee events (points), and (c) Z ee events (points), compared to Monte Carlo 
simulations (line), (d) Constraints on a and b from Z —> ee decays (large ellipse), J/ip —> ee 
decays (wide band), vr^ —> 77 ^ e'’'e“e"’“e“ decays (narrow band), and for all three combined 
(small ellipse). 


Z data and other data samples. Below is a discussion of the determination of the parameters 
in the W Monte Carlo. The Z data are treated in an analagous fashion. Systematic errors 
are set using large statistics samples of MC data and varying the parameter within its errors 
and are discussed throughout. 

The W production is modelled by the double differential cross section in priW) and 
rapidity, y, calculated at next-to-leading order by Ladinsky and Yuan [^] and using the 


MRSA pdf. The W resonance is generated by a relativistic Breit-Wigner, incorporating 
the mass dependence of the parton momentum distribution: 
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FIG. 5. (a) Comparison of the ti|| distribution from W ^ ev events (points) and the MC 

simulation (histogram); (b) Comparison of the angle between the recoil and the electron in the 
transverse plane from W ^ eu events (points) and the MC simulation (histogram). 


da e 

dm m (m2 - ' 


The angular decay products are generated at 0{as), allowing priW) > 0, in the W rest 


frame. This angular decay is of the form [14 


da 


dcosO 


(1 + ai cos 6^ + 012 cos^ 6) 


where Q!i ,2 = ai, 2 (PT(hT)) [0]- Radiative decays {W —*• eu'y) are generated according to 
Berends and Kleiss |^. Events in which W ^ tu ^ evvv are indistinguishable from 
IT —> ez/ decays and are therefore modelled in the simulation, including the polarization 
of the T in the decay angular distribution. The decay products are then boosted to the 
laboratory frame. At this point, the values of PtY) and Pt(IT) have been generated and 
is calculated. The effects of the detector and underlying event are now modelled. 

The EM (electron) calorimeter energy scale of the calorimeter was determined using 
J/Y —> ee, 7r° —77 —> e+e“e+e“, and Z —>• ee events. From test beam studies, it 
was determined that a linear relationship between the true and measured energies could 
be assumed: E^^as = tt-^true + This gives a relation Mmeas = OiMtme + 5f between 
the measured and true mass values, keeping terms to hrst order in 5 only. The variable 
/ = sin^ ^ depends on the event decay topology. Since the ratio of Mw to Mz is 

actually measured, one hnds 



(Mw'' 
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) 

1 + — 

{Mz — M\y) 

\Mz) 

V Mz . 

1 
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-1 

1 
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FIG. 6. (a) The transverse mass distribution, (b) the electron transverse momentum distri¬ 

bution, and (c) the neutrino transverse momentum distribution are shown for W events (points). 
The best fits of the MC simulation (histograms) are also shown. 

We note that to first order the measured ratio is insensitive to the EM energy scale, if 6 is 
small, and that the error on the measured ratio due to the uncertainty in 6 is suppressed. 

Figure ^ shows the mass spectra for the 7r°, J/il) and Z data samples. The allowed ranges 
for a and 5 are shown in Fig. |^d for each data sample. The overlap region is the la contour 
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FIG. 7. The D0 determination of the mass of the top quark is shown versus the measured W 
mass. The SM prediction (see text) for various assumptions of the Higgs boson mass is indicated 
by the bands. 


from all three data samples. The scale a is hxed by the Z data. The value of 6 is constrained 
by the J/ip and vr® data, essentially independent of a. Allowing a quadratic term in the 
energy response, to account for nonlinear responses at low energies as measured at the test 
beam, leads to the systematic error on 6. The allowed values determined for a and 5 are 
a = 0.95329 ± 0.00077 and 5 = —0.160 ± 0.016(sfaf.) ^; 2 io(s|/sf.) GeV. The error in the 
EM energy scale introduces an uncertainty in Mw of ±65 MeV/c^ and is dominated by the 
statistical error in determining the Z mass. 

The EM energy resolution is parameterized as a/E = \JC‘^ + {S/^JEj-Y ± [N/EY for the 
central calorimeter. Test beam data are used to set the sampling term, S = 0.13 (GeV^/^), 
and the noise term, N = 0.4 GeV. By constraining the width of the Z invariant mass 
distribution in the MG to that from the data, the constant term is set to G = (1.15 toJe). 
The uncertainty in the energy resolution leads to an uncertainty of ±23 MeV/c^ in Mw- 

The hadronic (recoil) energy scale of the calorimeter is determined relative to the EM 
energy scale by using Z events and measuring the transverse momenta of the Z from both 
the recoil or the two electrons. The balance is constructed: 


Pt— balance = [pT(ee) ±pT(rec)] ■ fj 


where fj is dehned as the bisector of the two electrons. From studies using HERWIG |T6 
and GEANT |^, it was determined that the recoil response could be written as a function 
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of the EM response: prirec) = TZrecPT{e.e) with IZrec = Orec + Prec^og priW). To ensure 
an equivalent event topology, Z events in which one electron is in the forward region are 
included in this study. Comparing data to MC in a plot of balance versus pr(ee) ■ p, the 
recoil response parameters are determined to be a^ec = 0.69 ± 0.06 and (3rec = 0.04 ± 0.02. 
The uncertainty in the recoil scale leads to an uncertainty of ±23 MeV/c^ in Mw- 

The recoil (hadronic) energy resolution is determined by modelling both components of 
the recoil to the W : pT{rec)^^°‘^ = IZrecPTi^^^c) + amb ■ Umb{tot) — 17(e). The hrst component 
is the “hard” component due to the initial pt of the boson. It is smeared using a Gaussian 
of width (Tree = Srec\JPTi^ec). The second component is the “soft” component due to 
the underlying event and is modelled by a minimum-bias event obtained from the data. In 
selecting the minimum-bias events to use, the luminosity distribution of the W event sample 
is modelled. The quantity Umb is the total of minimum-bias event and a^b is a scale 
factor. The amount of underlying event in the electron direction, 17(e), is subtracted from 
the recoil and added onto the electron momentum. Using the width of the balance 
distribution (to which the energy calibration has been applied), the values of Srec and amb 
are constrained. The measured values are s^ec = 0.49 ±0.14 and a^nb = 1.03 ±0.03 and their 
errors and lead to an uncertainty of of ±33 MeV/c^ in My/. 

Selection biases due to radiative decays and the amount of recoil energy in the electron 
direction and trigger efficiences are modelled in the MC simulation. The uncertainty in M/y 
due to the modelling of these efficiencies and biases is negligible in the £t to the Mt{W) 
spectrum. 

Backgrounds to the W event sample are included in the htting procedure by including 
the shape and fraction of background events. The largest source of background in the W 
sample is due to QCD multijet production in which there is a jet is mis-identihed as an 
electron and $j< due to energy fluctuations. This background contributes 1.4 ± 0.2% to the 
W sample. The other background considered is Z —ee events where one electron is not 
identihed. This background contributes 0.33 ± 0.06% to the W sample. The uncertainty in 
size and shape of the backgrounds gives an uncertainty in My/ of ±12 MeV/c^. All other 
sources of background are negligible. 

The last systematic error to consider is that due to the modelling of the W production. 
This uncertainty is due to the correlated uncertainties in the Pt{W) spectrum and the pdf’s. 
There are three phenomenological parameters in the production model calculation (pi, g 2 , Qs) 
m and the largest sensitivity of the pt spectrum is to the p 2 parameter. To constrain the 
production model, the gi and ps parameters are hxed to their nominal values and the value 
of p 2 is constrained by the Pt{Z) distribution from the data. Then the dependence of Mw 
on the pdf used in the theoretical calculation is measured from the difference in Mw from 
the nominal pdf (MRSA) as seen in Table For each pdf, the theoretical calculation uses 
the value of p 2 constrained by the data for that case. The uncertainties on the measured 
Mw due to the value of p 2 and the pdf used are ±5 MeV/c^ and ±21 MeV/c^, respectively. 
Errors on Mw are also ascribed to uncertainties in the value of Tve —> ± 9 MeV, the parton 
luminosity parameter —> ± 10 MeV, and the modelling of radiative decays —> ± 20 MeV. 
The total uncertainty on Mw due to the production model is a{Mw) = ±34 MeV. 

A measure of how accurately the MC describes the data is shown in Fig. The quantity 
u\\ = pT{rec) ■ e, which is the hadronic energy in electron direction, is shown in Fig. |^a. 
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pdf 

constrained §2 

AMw MeV 

MRSA 

0.59 

- 

MRSD- 

0.61 

+20 

CTEQ3M 

0.54 

+5 

CTEQ2M 

0.61 

-21 


TABLE IV. The variation of the measured W mass when using different pdf’s in the production 
model. Each theoretical calculation uses the constrained value of 52 for that pdf. 


A bias in u\\ directly affects the Mt{W) spectrum and it is also very sensitive to the recoil 
resolution. Another sensitive quantity is the difference in the azimuthal angle, A0, between 
the electron and the recoil and is shown in Fig. |b. Excellent agreement between the data 
and MC simulation is obtained. 


Source 

Statistical (W events) 

(Tm(w) in MeV/c^ 
69 

Statistical (Z events) 

65 

Non-Uniform energy response {rf) 

10 

Electron Angle Calibration 

28 

Electron Energy Resolution 

23 

Electron Energy Linearity 

20 

Electron Underlying Event 

16 

Hadronic Energy Scale 

20 

Hadronic Resolution 

33 

Pt{W) Spectrum 

5 

pdf 

21 

parton luminosity 

10 

W Width 

9 

Radiative Decays 

20 

QCD background 

11 

Z background 

5 

Systematic Total 

70 

Total 

118 


TABLE V. Summary of errors on the W mass measurement. 


The Mrp{W) distribution from the data is shown in Fig. |^a together with the distribution 
from the best fit value of Mw from the Monte Carlo simulation. The data are fit over a region 
60 to 90 GeV/c^ and the preliminary value of the W mass determined is Mw = 80.450 ± 
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0.070(stat.) ±0.065(sca/e) ±0.070(s|/st.) GeV/c^, giving a total error of ±118 MeV/c^. The 
errors on the W mass are detailed in Table 

As consistency checks, the prie) and Pt(i^) spectra are also fit to determine Mw as shown 
in Figs. ||b and |c. The fit to the pr(e) spectrum gives Mw = 80.49 ± 0.14 GeV/c^ and the 
Pt{i^) fit gives Mw = 80.42±0.18(stat.) GeV/c^ with the fitting region from 30 to 50 GeV/c 
in both cases. 

In summary, the measured W masses from the D0 data sample of IT —>■ eu decays 
with the e in the central r] region are = 80.35 ± 0.27 GeV/c^ (run lA) and = 
80.45±0.12 GeV/c^ (run IB, preliminary). Gombining these results and taking into account 
the correlated errors gives a D0 combined value of M\y = 80.44 ± 0.11 GeV/c^. 

Gombining the new D0 result with other measurements 


18,19 


from hadron collider 

experiments gives a new preliminary hadron collider average of Mw = 80.41 ±0.09 GeV/c^. 

The constraints placed on the Higgs mass can be seen in Fig. which shows the measured 
values of Mw versus rritop m compared to the Standard Model prediction for different 
values of the Mh in the rritop — Mw plane. 


V. CONCLUSIONS 

In conclusion, D0 has collected ~ lOOpb”^ of data from run 1 of the Tevatron and 
preliminary results on IF boson properties are found to be in agreement with the Standard 
Model. The IF and Z cross sections are measured in the e, /r, r decay modes. The IF 
width is measured to be Fvk = 2.159 ± 0.092 GeV. We confirm e — r universality in IF 
decays with the measurement /gY = 1.004 ± 0.032. The run 1 combined D0 IF mass, 
Mw = 80.44 ± 0.11 GeV/c^, is currently the most accurate direct measurement. 
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